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ABSTRACT 

We present a unifying empirical description of the structural and kinematic properties of all spheroids 
embedded in dark matter halos. We find that the stellar spheroidal components of galaxy clusters, 
which we call cluster spheroids (CSphs) and which are typically one hundred times the size of normal 
elliptical galaxies, lie on a "fundamental plane" as tight as that defined by ellipticals (rms in effective 
radius of ~ 0.07), but that has a different slope. The slope, as measured by the coefficient of the logo- 
term, declines significantly and systematically between the fundamental planes of ellipticals, brightest 
cluster galaxies (BCGs), and CSphs. We attribute this decline primarily to a continuous change in 
M e /L e , the mass-to-light ratio within the effective radius r e , with spheroid scale. The magnitude 
of the slope change requires that it arises principally from differences in the relative distributions of 
luminous and dark matter, rather than from stellar population differences such as in age and metal- 
licity. By expressing the M e /L e term as a function of a in the simple derivation of the fundamental 
plane and requiring the behavior of that term to mimic the observed nonlinear relationship between 
log M e / L e and logc, we simultaneously fit a 2-D manifold to the measured properties of dwarf ellipti- 
cals, ellipticals, BCGs, and CSphs. The combined data have an rms scatter in logr e of 0.114 (0.099 for 
the combination of Es, BCGs, and CSphs), which is modestly larger than each fundamental plane has 
alone, but which includes the scatter introduced by merging different studies done in different filters 
by different investigators. This "fundamental manifold" fits the structural and kinematic properties 
of spheroids that span a factor of 100 in a and 1000 in r e . While our mathematical form is neither 
unique nor derived from physical principles, the tightness of the fit leaves little room for improvement 
by other unification schemes over the range of observed spheroids. 

Subject headings: galaxies: formation — galaxies: elliptical and lenticular, cD — galaxies: fundamental 
parameters — galaxies: structure 



1. INTRODUCTION 

Given the differences in how, where, and when galaxies 
form, there are some relationships among galaxy param- 
eters that appear impossibly tight. The importance of 
such relationships is that they must evince deeper and 
yet unknown laws of galaxy formation. The principal 
example of such a relationship is that relating scale r e , 
surface bright ness I e ., and velocity dispersion a for ellip- 
tical galaxies ({Dressier et al.l Il987t iDiorgovski fc Davisl 
I1987D . In the 3-space defined by these parameters, el- 
liptical galaxies lie on a plan e with modest scatter (~ 
0.1 dex: iJorgensen et alllToM iBernardi et al.ll200.^ . Al- 
though the origin of this remarkable result is not fully un- 
derstood, this relationship is widely known as the "Fun- 
damental Plane" (hereafter, FP). 

There have been previous efforts to place all dynam- 
ically hot systems onto the FP, or related r elationships 
(e.g., lKormendvlll985l iBurstein et"aTl I1997H . Althou gh 
systems from gl obular clusters fl Burstcin et al.lll997l) to 
galaxy clusters l)Schaeffer et al.lll993ji do obey general 
relationships similar to the elliptical FP, the slopes and 
intercepts are often significantly different, suggesting dif- 
ferences in the role of dissipation, the distribution of 
angular momentum, and the relative distributions of 
light and dark matter. Although we have yet to un- 
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tangle the clues provided by these relationships, probing 
the extremes of the family of FPs — where spheroids 
must arise from qualitatively different formation pro- 
cesses — may provide a breakthrough. For example, 
while dissipation is clearly important in boosting the 
central phase space density of normal Es l)LauerlH985t 
lCarlbergi n"9861 . the hierarchical accretion that produces 
the largest spheroids, the intracluster ligh t component of 
galaxy clusters, is mostly dissipationless llMurante et alJ 
l2004t IWillman et alJl200l iSommer-Larsen et alJ l2005). 
Differences may then arise in how spheroids across these 

scales populate (r P _, J e , q)-s pace. 

Originally identified bv iMatthews et alJ l|1964fl . the 
light in excess of the standard r 1 / 4 surface brightness pro- 
file observed at large radii in some brightest cluster galax- 
ies (BCGs) is generally referred to as the "cD envelope" , 
although specialists have a wider range of term inology 
(IOemlerlll97l fl97fit ISchombertlll98fi psl llTlsl . More 



recent work llUson et all 



19M l199lVlScheickfcKuhn 



et 

119941 iGonzalez et al.l 12000 "1 questioned the existence of 
such a component and suggested that difficulties in back- 
ground subtraction at such low surface brightness lev- 
els might h ave skewed measurements. However, the lat- 
est e fforts llFeldmeier et all 120021 120041 IGonzalez et~aT] 
l2T)0llabet"ti et alJEOOSlT have all found excess luminos- 
ity above the r 1 / 4 profile of BCGs, although at a surface 
brightness well below that associated with classical "cD 
envelopes" . This luminosity is attributed to intraclus- 
ter stars (hereafter, ICS), a stellar population that is 
gravitationally bound to the cluster but unbound to any 
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individual cluster galaxy. The ICS component is well de- 
scribed by an r 1 / 4 surface brightness profile and satisfies 
a relationship between r e an d J R that is similar t o that 
of normal elliptical galaxies ( Gonzalez ct al. 2005, here- 
after Paper I) . As such, the stellar component associated 
with the cluster as a whole appears to be an extremely 
large version of normal spheroids, so we designate it as 
the cluster spheroid (CSph). 

To examine the place of CSphs among more traditional 
spheroids requires a measurement of the CSph velocity 
dispersion in clusters where we have already observed 
the ICS (Paper I) . Although the ideal experiment would 
involve measuring the velocity dispersion of the ICS di- 
rectly, that tas k is quite difficult be cause of its low surface 
brightness (see lKelson et al.ll2002|) and is not possible for 
our large sample of clusters without a tremendous invest- 
ment of large telescope time. Therefore, we adopt the 
cluster galaxy velocity dispersion as the CSph dispersion 
and discuss possible shortcomings of this approach. We 
review briefly the status of the ICS measurements and 
present our cluster kinematic measurements in I n §3 
we discuss the connection between the BCG and ICS, 
the FP populated by the CSphs, and a new relationship 
between the FPs of ellipticals, BCGs, and CSphs. In §5 
we summarize our findings. 

2. DATA 

In Paper I we present observations of 24 clusters and 
rich groups that span a range of velocity disp ersions and 
Bautz- Morgan types l)Bautz fc Morganlll970j) . The sam- 
ple consists of nearby systems (0.03 < z < 0.13) that 
contain a dominant BCG with a major axis position an- 
gle that lies within 45° of the east-west axis (the drift 
scan direction). We present details of these unique data 
and the reduction procedure in Paper I. Here we discuss 
the 23 of those clusters for which we obtained kinematic 
data (Abell 2376 is the one cluster from Paper I that 
is omitted here). All calculated quantities in this pa- 
per are for a concordance cosmology model: fl m = 0.3, 
fl A = 0.7, and H = 70 km sec" 1 Mpc -1 . 

We obtained spectra of galaxies in the fields of the 12 
clusters from Paper I with insufficient kinematic data 
in the literature. We observed those galaxies with the 
Hydra spectrograph at the CTIO Blanco 4m telescope 
on the nights of 28 July 2003 to 31 July 2003 using the 
KPGL1 grating, which has 632 lines mm -1 and provides 
a dispersion of 0.59 A pixel -1 . The measured spectra 
span 3700 to 6084 A. In combination with the bare fibers 
(no slit plate was used), the resulting spectral resolution 
is 4.1 A. Targets were selected from our optical imaging 
(described in Paper I) to satisfy the simple criteria that 
they be fainter than the BCG and lie within a projected 
distance of 1.5 Mpc from the BCG. Priority for fiber as- 
signments was based on galaxy magnitude. We obtained 
flat fields and wavelength calibrations through the fibers 
at each target position and typically observed each clus- 
ter field for 1800 sec, three to four times. Total exposure 
times vary due to variable weather conditions, but the 
above numbers are representative. 

The data were reduced in a standard manner, includ- 
ing bias subtraction using both overscan regions and bias 
frames, and flat fielding using continuum lamp exposures 
taken through the fibers. We used these same flat field 
exposures to define and trace the apertures, allowing a 



spatial shift on the detector to best fit the position of 
the spectra in each individual exposure. We measured 
redshifts using cross-correlation techniques, specifically 
the IRAF 4 tasks within the RVSAO package, XCSAO 
and EMS AO. We visually inspected all spectra both to 
reject marginal cases and to add systems with emission 
line redshifts that were missed by the cross-correlation. 
Most absorption line redshifts come from cross correla- 
tions with R > 4. For spectra with both absorption and 
emission line redshifts, we give the average of the two 
measurements. For the remaining 11 cluster fields from 
Paper I that were not targeted spectroscopically, we ob- 
tained redshifts from the literature (using the NASA Ex- 
tragalactic Database, NED) for galaxies projected within 
1.5 Mpc of the BCG. 

The kinematic results from our observations and lit- 
erature search are given in Tabled To determine clus- 
ter mem bership, we employ a pessimistic 3er clipping al- 
gorithm llYahjlfe Vidallll977ft using biweight estimators 
l)Beers et alJ 119901) of location and scale for the cluster 
mean velocity and velocity d ispersions, respectively (cf. 
iZabludoff fc Mulchaevi ri998). The resulting number of 
cluster members, their velocity range, mean velocity, and 
line-of-sight velocity dispersion are presented in Tabled 
with the source entry referring to whether the data come 
from our own Hydra (H) observations or from NED (N). 
Throughout the rest of this study, we use the line-of- 
sight velocity dispersion for the velocity dispersion of the 
CSphs, assuming that the conversion factor is the same 
for all clusters and hence absorbed into the constant term 
in log — log relations (see §3.2, 3.4, and 3.5). 

Finally, we present the total luminosity of the resolved 
galaxies within the ICS r e in Table ^ We include this 
measurement because we will explore whether various 
relationships are strengthened when using the entire lu- 
minous content of clusters. We compute a background- 
subtracted luminosity within r e , including all galaxies 
fainter than th e BCG down to m = 19.5. We use 
the SExtractor ijBertin fc Arnoutsil99 6) stellarity index, 
which is robust down to this magnitude, to exclude stars, 
and compute the background level using galaxies located 
at distances of 23.3' - 46.6'(2000-4000 pixels) from the 
BCG. We also apply a completeness correction to ac- 
count for the contribution of faint cluster galaxies at 
m > 19.5. For this correction we adopt a = — 1.21 
and AU.r = -21.96 (jChristlein k Zabludofl1l2003|) and 
R — I = 0.82 to convert from R to I. 

3. RESULTS AND DISCUSSION 

3.1. Effective Radius and Surface Brightness 

A key to the existence of the FP is the inverse rela- 
tionship between r e and I e , the effective radius and the 
mean surface brightness within that radius, respectively. 
In Paper I, we demonstrate that both the inner compo- 
nent of brightest cluster galaxies, which we associate with 
the BCG itself, and the outer component, which we refer 
to as the intracluster light or intracluster stars (ICS), ex- 
hibit a tight relationship between r e and I e , although the 
slopes of these relationships for the BCG and ICS com- 
ponents differ slightly. We adopt the general practice of 
expressing r e in units of kpc and I e in units of Lq/pc 2 . 

4 IRAF is distributed by the National Optical Observatories, 
which are operated by AURA Inc., under contract to the NSF. 



Fundamental Manifold of Spheroids 



3 



TABLE 1 

Measured Cluster Parameters 



Cluster 


M IC ,s 


r e,ICS 

[kpc] 


M SC G 


r e,BCG 

[kpc] 


M Bi GAL 


Source 


No. 


Range 


V 

[km s- 1 ] 


<j 


A0122 


-25.56 


107.9 


-22.97 


4.7 


-24.97 


H 


28 


[32410,35611] 


34020 ± 130 


680+"° 

— 90 




±0 01° 9£ L 


±2 9 9 '2 , 


±0 03°'°5 

xu ' UJ -0 .05 


±0 I 2 - 
J -"' 1 - 0.2 












A1651 


-24.95 


477.2 


-25.21 


48.2 


-25.76 


N 


28 


[23338,27324] 


25400 ± 190 


990+" 9 

— 100 




±0 07°- 4S „ 


±53 I 23 , 9 ,! . 


±0 01° S°„o 

_1_U.UX_Q Q2 


±0 5°-£ 

_i_u. j_0 2 












A2400 


-24.61 


184.6 


-24.19 


20.4 


-24.37 


H 


50 


[25147,27660] 


26380 ± 90 


650+°° 

— 60 




±0 04°-i 9 

xu,w - 0.44 


±44 9 169 ; 7 
- 1 - ^— 75.2 


±0 08°'J 7 1n 
- LU - wo — 0.11 


±1 0° ? j 
xi,u - 1.4 












A2401 


-24.48 


27.6 


-22.30 


1.8 


-20.13 


N 


24 


[16123,18207] 


17120 ± 100 


460+1°° 

— oU 






±0 4 2 ? „ 

- Lu- ^— 1.9 


±0 03°i°.., 


±0 I 2 

-LVJ. -L_ o.2 












A2405 


-23.38 


93.9 


-22.89 


3.8 


-19.09 


H 


8 


[10778,11308] 


11020 ±70 


140+|° 

— dO 




-^■^ — 0.40 


±4 7 6? o? r 

' -33.5 


±0 01°'S 6 n B 

_1_U.U-L_Q Qg 


±0 I 2 












A2571 


-25.11 


107.7 


-23.89 


5.4 


-21.49 


H 


11 


[31103,34517] 


32710 ± 110 


670 + ?n 
— 70 




±0 01°'J 7 ^ 

IU ' ul -0.22 


±5 43 n5 n 


±0 02° S S „„ 


±0 I 4 „ 
=EU * X — 0.4 












A2721 


-25.02 


47.9 


-23.14 


18.1 




N 


67 


[32452,36859] 


34300 ± 100 


840+8° 

— 70 




±0 03° S 4 n, 

- LVJ,WlJ — 0.04 


±0 8 1 -? , 


±0 14 - 08 
-^■^ — 0.08 


±1 6 1 ? , 
xl,u - 1.1 












A2730 


-26.00 


245.9 


-23.85 


8.1 


-25.18 


N 


13 


[34281,37767] 


36000 ± 300 


1020 + l^n 
— 140 




±0 03°-R 5 ,^ 


±13 5 12 , 1 ,; 9 , 


±0 03° i 9 1c 


±0 2 1 ? , 
- Lu- ^— 1.3 












A2811 


-25.53 


139.5 


-23.59 


5.5 


-24.72 


N 


35 


[30365,34227] 


32350 ± 150 


860+"!i 

— 100 




±0 03°- 22 „„ 


±8 6 59 ,I . 

— 44.4 


±0 03 22 .. 


±0 2°- 8 „ 












A2955 


-25.22 


93.5 


-22.61 


3.4 


-24.24 


H 


22 


[27539,29118] 


28260 ± 70 


320+ 9 ° 

— 70 




±0 oi - 08 ™ 


i2.3_g_g 


±0 03°'S 9 n 


±0 , 

J - u ' 1 - 0.3 












A2969 


-25.35 


232.0 


-23.90 


19.7 


-24.64 


N 


21 


[35684,40268] 


37740 ± 220 


980+{ 9 ° 

— loO 




±0 03°- 27 .„ 


±26 7 20 o 3 o 6 o 

_i__,u. I _g 3 2 


±0 lO -;. 1 ,,. 


±0.8 1 f 7 












A2984 


-25.79 


202.7 


-23.04 


6.1 


-25.22 


H 


29 


[29882,32379] 


31230 ± 90 


490+"° 

— 90 




±0 Ol^S 1 ™ 


±5 4 23 o 2 „ 

— 2.4 


±0 02° pi 


±0 I 2 „ 

_l_U. -L_o 












A3112 


-25.74 


102.1 


-22.95 


8.3 


-23.79 


N 


59 


[19957,24560] 


22560 ± 120 


940+T 4 ° 

— 120 




±0 Ol 06 


±1 7 12 - 4 


+0 04 16 


±0 2°- 9 












A3166 


—24.65 


69.2 


—21.92 


1.7 


—21.73 


H 


12 


[34794,35511] 


35190 ± 70 


210+ 40 
40 




+0 02 07 




in n4 0.08 
3ZU.U*_q g 


+0 I - 2 












A3693 


—25.46 


314.5 


—23.56 


4.9 


—24.81 


H 


31 


[33500,38657] 


36840 ± 190 


1030+™ 
XUJU — 130 




_i_n (150.20 


+24 46L5 


+0 02° 10 


+0 1°- 3 












A3705 


—24.27 


30.8 


—21.99 


2.1 


—22.57 


N 


53 


[25110,29106] 


26820 ± 140 


1010+?° 




+0 m 003 


+1 o 10 - 9 


_l_n n7 0.02 
=tU.U. _ 25 


+0 1 15 












A3727 


-24.97 


215.2 


—23.07 


6.1 


—24.16 


H 


25 


[33295,36030] 


34760 ± 120 


580_ 1 qo 




_i_ n n^O iS 


+17 1 54 - 6 


+0 mO.os 


+0 1°- 2 












A3809 


—23.85 


126.8 


—24.14 


20.8 


—22.79 


N 


69 


[17550,20192] 


18680 ± 70 


540+ 60 




±0 03°i 3 ,,. 


±17 2 31 4„ 
__:._ 1 -^_22.0 


±0 03° R 2 „^ 
xu,ud - 0.02 


±0 4 01 , n 
:cu '^— 0.2 










A3920 


-24.80 


45.6 


-23.48 


5.4 




H 


17 


[37120,39232] 


38000 ± 110 


430+ 14 ° 

— 1UU 




in mll.Ol 
ItU.Ul_ 04 


±3.2 4 _ 4 2 


+0 10 07 
_tu.±u_ 05 


±0-4^.2 












A4010 


-25.56 


186.3 


-23.52 


13.7 


-24.49 


N 


31 


[26957,30069] 


28620 ± 120 


630lf 2 ° 




±0.02°_i 9 2g 


±10.7^ 7 g 4 


±0.06°_ 37 2g 


±0.6 2 4. 












APMC020 


-24.34 


180.2 


-24.40 


17.6 


-23.41 


H 


10 


[32477,33228] 


32940 ± 80 


240l«° 




±0.040_-25 62 


±34.4 2 f i i 3 




±0.4°_- 9 7 












AS0084 


-25.33 


72.0 


-21.94 


2.0 


-23.60 


N 


24 


[31436,34068] 


32390 ± 120 


520+ 160 




_l_n niO.03 
=tU.Ul_ 05 


+ 1 6 5 - 4 


__0.04°_° 8 10 


+0 1 - 2 












AS0296 


-24.89 


46.1 


-23.99 


10.0 


-20.36 


H 


34 


[20082,21872] 


20860 ± 90 


480+^ 




+0 01 °- 01 


±l-2 7 -1.x 


±0.04°_ 39 12 


±0.3°_- 9 3 













Because the luminosity within the ICS r e comes from 
at least three components (BCG, ICS, and resolved 
galaxies), there is some ambiguity in deciding which lu- 
minosity to use in the calculation of r e and I e for the 
CSph. In all four panels of Figure ^ we choose to plot 
r e for the ICS component, but alter our choice of the 
luminosity used to calculate I e . In the upper left panel, 
we include only the surface brightness of the ICS. By 
construction for a r 1 / 4 profile, half of the total ICS lumi- 



nosity is enclosed within r e . In the upper right panel, we 
have added the entire BCG luminosity to that of the ICS 
enclosed within r e . The addition of the two luminosities 
means that the adopted r e no longer encloses half the 
total light. However, the BCG component is typically 
a small fraction of the total light of the combined sys- 
tem (Paper I) and so does not drastically affect r e , i.e., 
the half-light radius for the combined two components is 
typically 20 to 30% smaller than the ICS r e . In the few 
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1.5 2 2.5 1.5 2 2.5 

log r e [kpc] 



Fig. 1. — The relationship between r e and I e for different def- 
initions of I e . In all panels we use r e measured from the ICS. In 
the upper left panel we use only the ICS luminosity within r e to 
calculate I e . In the upper right panel, we add the BCG luminosity 
to the calculation of I e . In the lower left panel, we add both the 
BCG luminosity and that of the resolved cluster galaxies within 
r e . In the lower right panel, we only use the luminosity of resolved 
cluster galaxies within r e . 



Fig. 2. — The relationship between r e and I e for both ICSs and 
BCGs. In both panels the open circles represent the data for the 
BCGs. In the left panel, the filled circles represent the data for the 
ICS. In the right panel, the filled circles represent the calculated 
I e using the sum of the ICS and BCG luminosities within the ICS 
r e . The solid line illustrates the direction along which errors in the 
measurement of r e will scatter the data. 



cases where the BCG component dominates, the com- 
bined light r e is several times smaller than the ICS r e . 
Because 1) the ICS generally dominates, 2) r e and I e are 
correlated, and 3) this correlation tends to slide the data 
along the r e — I e relationship when r e changes, our use 
of the combined light r e , instead of the ICS r e , does not 
significantly alter the slope or scatter of the various rela- 
tionships we discuss. Regardless of which r e is adopted, 
including the BCG component in the calculation of I e 
decreases the scatter. Lastly, we include the luminos- 
ity from resolved cluster members within r e (lower left 
panel of Figure The empirical result, for which we 
will discuss a possible interpretation below, is that the 
r e — I e relation is tightest when we combine the luminos- 
ity from the central BCG and the ICS to describe the 
CSph. This result is supported by a quantitative mea- 
sure of the scatter and does not rest on a few deviant 
points. When the BCG luminosity is included, the eight 
clusters with a residual that is > 0.1 in the ICS r e — I e 
correlation all drop to below 0.1 and no clusters increase 
their residual to > 0.1. 

A difficulty in interpreting the r e — I e relationship is 
that I e depends directly on the determination of r e . This 
degeneracy does not render the relationship meaningless 
because there is additional information in I e regarding 
the luminosity profile. However, the degeneracy does im- 
ply that errors in r e will translate into correlated errors 
in I e . The observed relationship is close in slope to that 
arising from correlated errors and is therefore somewhat 
suspect (Figure . At this point one must either have 
confidence in the quoted errors or provide independent 
measurements demonstrating that r e and I e are not spu- 
rious. Although our uncertainties are indeed sufficiently 
small that we claim not to have confused a large r e ICS 
component with a small one (see the original Figure 9 



in Paper I or Table QJ, these are difficult measurements 
that are susceptible to systematic errors, particularly in 
the sky determination (see Paper I). Given both the ran- 
dom and systematic errors, we cannot easily rule out the 
possibility that some luminosity is artificially taken from 
one component and placed into the other by our fitting. 

We first examine whether another measurement can 
indirectly validate our measurement of r e . For example, 
our measurement of r e for the ICS correlates with the 
velocity dispersion of the cluster, a, sufficiently well that 
we can rule out at the 97% confidence level that such 
a correlation arose randomly. This correlation suggests 
that there is real information in the r e measurements (a 
lack of a correlation would not conversely suggest that 
the r e values were dominated by errors). In contrast to 
this ICS result, the r e values for the BCG correlate only 
weakly with a and there is a 13% chance that such a 
correlation could arise at random. 

Next we examine whether the measurement of the ICS 
luminosity correlates with cluster velocity dispersion. We 
find again that there is a strong correlation (98% confi- 
dence level). Given the large uncertainties in a (Fig- 
ure [3J , the underlying correlation must be significantly 
tighter, which once again suggests that errors in the sky 
determination are not sufficiently large to disrupt our 
measurements of the ICS. As with r e , the correlation for 
the BCG component is much weaker (a 20% chance of 
random) . 

Finally, we compare the ICS magnitude to that of the 
summed cluster galaxies inside r e (Figure . We find 
a strong correlation (99.97% confidence level) between 
these two quantities, and again the correlation is much 
weaker between the BCG and the resolved galaxy lumi- 
nosity (27% chance of random). We conclude that our 
measurements of r e and L e , which in turn define I e , for 
the ICS are reliable and that the range of these values 
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Fig. 3. — Cluster velocity dispersion versus ICS magnitude. Er- 
rorbars in Mjqs represent the potential systematic error due to 
a lo" systematic sky error. The formal fitting errors in Mjqs are 
significantly smaller than the plotted systematic errors. 
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Fig. 4. — Magnitude of summed cluster galaxies inside r e versus 
the magnitude of the ICS. 



does have physical meaning. Nevertheless, when deriving 
quantitative measurements we take care to account for 
the correlated errors through Monte-Carlo simulations. 

We conclude from these tests that our structural mea- 
surements, particularly those for the ICS, accurately re- 
flect the properties of the cluster. We are left with two 
families of interpretations for the measured improvement 
in the r e — I e relation when the combination of the BCG 
and ICS luminosity is used. The improvement may sim- 
ply arise because our partition of luminosity between the 
BCG and ICS is poor, and therefore the correlation using 
the sum of components has less scatter than that using 
one poorly determined component. Alternatively, the 
improvement may indicate some deeper connection be- 
tween the BCG and ICS and that the two together prop- 
erly reflect the CSph. We favor the latter because (1) 



adding the BCG luminosity reduces the scatter both for 
systems with too little and too much luminosity within 
r e , (2) there is no relationship between the scatter in 
the r e — I e relation and the estimated systematic error, 
and (3) the decomposition of components in Paper I was 
based not solely on a surface brightness profile, but also 
on ellipticity and position angle variations. We also favor 
the latter because there are external arguments for a con- 
nection between the properties of BCGs and the global 
cluster environment. These connections include the find- 
ing that BCGs are not drawn fro m the luminosity func- 
tion of other cluste r ellipticals l)Tremaine fc Richstonel 
H977tlDresslerlll97^ . and that there is a correspondence 
between the BC G optical luminosity and the cluster 
X-ray luminosity llHudson fc Ebelindll99''^: iBurke et alJ 
120001 INelson et al.ll2002() . In the final analysis, none of 
these arguments is definitive, and one could prefer to sim- 
ply fall back to the empirical improvement in the scatter 
to adopt the joint BCG and ICS luminosity. While we 
choose to use that combination as the best description 
of the CSph, we also present results using only the ICS 
luminosity, and show that this choice does not materially 
affect any of our findings and conclusions. Unless other- 
wise noted, from now on we use the sum of the BCG and 
ICS components for the CSph. 

3.2. Defining the CSph Fundamental Plane 

Emboldened by the strong r e — I e relationship for the 
CSphs, and following the example of the FP for elliptical 
galaxies, we examine whether the cluster velocity disper- 
sion helps further refine the portion of parameter space 
occupied by our cluster spheroids. While the most ap- 
propriate analog would be the velocity dispersion of the 
ICS component itself, the low surface brightness of the 
ICS precludes measurem ent of a in all but a few clusters 
fsee iKelson et all i|2002j) '). We showed above that vari- 
ous properties of the ICS are correlated with the cluster 
velocity dispersion as measured from the galaxies, and 
so the use of the cluster velocity dispersion as a proxy 
for the dispersion of the ICS is not without reason. The 
data are so far inconclusive on the question of whether 
this is a good assumption. In one cluster the ICS dis- 
persion is consi stent with that me asured from galaxies 
at la rge radius IKelson et al.ll2002|) . and in another it is 
not l)Gerhard et al] 120051) . A further departure of out- 
work from a direct analogy with the study of ellipticals 
is that the velocity dispersions of ellipticals are measured 
within r e , or corrected to represent the dispersion within 
r e . Because there are typically few cluster galaxies within 
r e , we adopt the global cluster velocity dispersion as our 
velocity measurement with no correction. All of these 
departures should weaken any underlying CSph FP rela- 
tion, and as such our results will be conservative. 

Before proceeding to the FP analysis, we ask if a is 
correlated with the residuals of any of the correlations 
examined so far. If so, such a correlation would further 
suggest that the cluster velocity dispersion is related to 
the structural properties of the CSph and provide addi- 
tional motivation to proceed with an examination of the 
FP. In Figure we show the r e — L e relationship, and 
the residuals about that relationship, A, versus a. We 
examine residuals about the r e — L e relationship rather 
than about the r e — J e relationship because r e and L e 
(= L/2) are independently measured quantities. The 
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Fig. 5. — Connecting r e , L, and a for the CSphs. Left panel 
shows the correlation between ICS r e and CSph (ICS+BCG) mag- 
nitude. The solid line represents the best fit line to those data. 
The right panel shows the correlation between the cluster veloc- 
ity dispersion a and the residuals to the best fit line from the left 
panel. 



Fig. 6. — Fundamental Plane for CSphs. We show projections 
of the FP for the ICS only (left panel), ICS+BCG (middle panel), 
and ICS+BCG+Galaxies (right panel) components. The line is 
the 1:1 line, which is the plane in these appropriately scaled axes. 
The numerical values for A, B, and C are given in Table l2l 



A — a relationship is highly significant (99.2% confidence 
level) and demonstrates that the velocity dispersion does 
provide additional information on the structure of these 
systems. 

The strong relation between the ICS effective radius 
and the combined mean surface brightness of the CSph 
(upper right panel of Figure ^ with a scatter of 0.15 
in the plotted units) already suggests that the principal 
extension of the CSph fundamental plane, if there is one, 
is primarily along these two axes. As is generally done 
for ellipticals, and to facilitate comparison with them, we 
include a as the third axis and fit 

AlogCT + Slog/ e -logr e + C = 0. (1) 

The fitting of a plane to these data is complicated 
for several reasons. First, in a typical fitting problem, 
the independent variables are much less uncertain than 
the dependent variables, and hence errors in the inde- 
pendent variables can be safely ignored. In this case, 
however, not only are the uncertainties comparable for 
all three variables, but the uncertainties in two of these 
variables are highly correlated. Second, there are multi- 
ple ways in which one can define residuals about a plane. 
In the direct method (for a more detailed d iscussion of 
the nomenclature see Bernardi et alJ l{200 3)). one min- 
imizes the residuals along an axis (typically, residuals 
in r e are quoted), while in the orthogonal method one 
minimizes the 3-D distance between the data and the 
plane. These two methods appear to give roughly the 
same answer for B in Eq. 1 when fitting Es, but dif- 
fer by about 25% in their determination of A in Eq. 1 
l)Bernardi et al .1120031) . Although the choice of method is 
important for quantitative measurements, the prime con- 
sideration when comparing results from various studies is 
whether the same method is used throughout. We choose 
to focus on results from the direct method for continuity 
throughout the paper and for the technical reasons that 



develop below. Therefore, we will refit some previously 
published data that were fit originally with the orthogo- 
nal method. For our exploration of the CSph FP alone, 
we present both direct and orthogonal fits in Table [21 

Although the fitting proble m has been solved an alyt- 
ically for Gaussian errors fsee iBernardi et aT]l2003j) . we 
opt for a numerical fitting method so that we can ex- 
periment with asymmetric error distribution, and even- 
tually move beyond planar fits in a straightforward way 
(§3.4). We fit a plane to our data, using r e as measured 
from the ICS surface brightness profile, I e as calculated 
from the adopted luminosity profile (ICS, ICS+BCG, or 
ICS+BCG+galaxies, see below), and a as measured from 
member galaxies. To account for measurement uncer- 
tainties, and the correlated errors of r e and I e , we ran- 
domly draw errors for r e , L, and a from the uncertainties 
given in Table Q] an d then recalculate L e and I e . We use 
the estimated systematic uncertainties resulting from a 
la global error in the background sky determination be- 
cause these are larger than the statistical errors for the 
measured parameters. Where errors are asymmetric, we 
adopt Gaussians of different dispersions in the two direc- 
tions. We repeat the fitting 1000 times and then deter- 
mine uncertainties in the fitted parameters by identifying 
the range that encloses two-thirds of the results from the 
1000 trials. These fits, particularly in cases where the 
errors are asymmetric, will not result in the same fitted 
parameters as a direct fit to the original data because 
the mean value of the distribution is not equal to the ob- 
served value. In general, the resulting difference between 
the fitted parameters is undetectable. 

In Figure El we plot edge-on projections of the fitted 
plane and data for the various plausible choices of the 
CSph luminosity. Again, as in §3.1, we explore adopting 
either the luminosity from just the ICS, from the combi- 
nation of the BCG and ICS, or from the combination of 
the BCG, ICS, and the resolved cluster members within 
r e . The scatter about the plane is 0.106 when we use 
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TABLE 2 

CSph Fundamental Plane Coefficients 



Method Components A B C r e RMS 



Direct 



Orthogonal 



ICS 0-26±nog -0A8±°°l 1.661^8 al06 

ICS+BCG 0- 21 iao2 -0.56±° Q ° \ 1.96±g;g| 0.074 

ICS+BCG+Galaxies 0.32± ;°£ -0.65l°;°| 1.86±°;^ 0.101 



ICS 0A7t° ° 7 7 — 0.51±gg| 1.10±°'?£ 0.112 



ICS+BCG 0.23I";"J -0.57IS'oi L91 -o'.g? °- 073 

ICS+BCG+Galaxies 0A9±° n °l -0.7l±° n °3 lAA±° n ^° 0.107 



only the ICS luminosity to determine I e , 0.074 when we 
use the BCG+ICS luminosity, and 0.101 when we also 
include the resolved galaxies within r e . Once again the 
combined BCG+ICS luminosity provides the tightest re- 
lationship 5 . Therefore, we conclude that the CSph is best 
described by the BCG+ICS and adopt this description 
for the CSph unless otherwise noted. None of the basic 
results we present depends on which definition is used 
(see Table0). The change in the scatter between the 0.15 
measured for the r e — J e relationship and the 0.074 mea- 
sured for the CSph FP demonstrates the improvement 
introduced by utilizing the velocity dispersion. The low 
scatter relative to the FPs of other spheroids suggests 
that using the velocity dispersion of the ICS rather than 
that of the galaxies is unlikely to significantly improve 
the relation. This result does not necessarily imply that 
the velocity dispersion of the ICS is the same as that of 
the cluster galaxies, but rather that they scale propor- 
tionally so that the scaling between the two is absorbed 
into the constant term in the FP equation. The fitted 
parameters and rms scatter between the fitted and ob- 
served log r e values are listed in Table |21 

3.3. Comparison to Other Spheroids 

The existence of a FP for CSphs leads naturally to 
the question of whether this FP is similar to that found 
for other, smaller spheroids. Due to a wide range of 
studies, data exist for systems as diverse as dSphs, dEs, 
Es, and BCGs. Of these, the most extensive studies of 
the FP have been done with normal Es, so we begin with 
a comparison to those galaxies, then compare to other 

5 We confirm that the reduction in the scatter when the BCG 
light is included is statistically significant by fitting 1000 random- 
ized versions of the data. We draw the relevant parameters from 
Gaussian distributions according to their estimated uncertainties. 
In this manner we are able to answer the questions: 1) if the tight 
correlation seen in the ICS+BCG data is correct, what is the like- 
lihood of getting as "poor" a correlation as that seen for the ICS 
data alone by random chance, and 2) if the noisier correlation seen 
in the ICS data alone is correct, what is the likelihood of getting 
as "tight" a correlation as that seen for the ICS+BCG data by 
random chance. From these Monte-Carlo simulations, we find that 
none of the 1000 trials of the ICS+BCG data produces as large a 
X 2 value as that of the fit to the ICS data alone and that none of 
the 1000 trials of the ICS data alone produces as small a \ 2 as that 
measured for the ICS+BCG data. 



spheroids, and finally attempt to place all spheroids on 
a single relationship. 

The degree of scatter in the CSph FP is slightly 
lower t han that typical l y obt ai ned for elliptical gala x- 
ies (see Uorgensen et all l|1996j) : iBernardi et alJ l|2003ft ). 
To enable a more direct comparison, we apply our di- 
rect fitting technique to the elliptical data compiled by 
Dorgensen et all (^)96). The direct fitted parameters 
(A = 1.00, B = -0.75, and a scatter of 094) d iffer 
slightly from those derived bv I.Torgensen et al |)199(tD us- 
ing the orthogonal fitting method (A = 1.24, B = -0.82 
and a scatter of 0.084), but we recover those published 
values when we fit the data using the orthogonal ap- 
proach. 

Both the direct and orthogonal fits to the elliptical 
galaxy data differ significantly from the corresponding 
fits to the CSphs, and that difference is primarily in the 
value of A. Therefore, while the CSphs fall onto a FP 
as tightly as do elliptical galaxies, it is not the same 
plane. This comparison, however, is complicated by the 
use of data obtained with di fferent filters (for example, 
r from lJor gensen et al| l|199 6f) an d Ic from Paper I). We 
correct the lJorgensen et all l|1996|) data to Ig, hereafter I, 
using a model color calculation for Es bv lFukugita et alJ 
(1995), but there remains some discrepancy between the 
two samples that depen ds on the degree to which (1) r e 
is color dependent (see iPahre et alJ l)1998l) for evidence 
of color dependence) and (2) the colors of all ellipticals 
differ from the models. Nevertheless, these uncertainties 
will only increase the scatter in any attempt to place the 
full range of spheroidals on one relati onship. 

We fit a plane to the BCG data of lOegerle fc Hoessell 
(1991), again correcting using the colors of e l liptica l 
galaxies and the calculations bv iFukugita et all |l995). 
In this case, we find a planar fit that has an A value 
of 0.52, which is between those of the ellipticals (1.00 
or 1.21) and CSphs (0.21). The BCG scatter about the 
best-fit plane (0.085) is comparable to that found for 
both Es and CSphs. These results suggest that each of 
these distinct types of spheroids (Es, BCGs, and CSphs) 
lies on a comparably thin, yet distinct, plane in this 
space. 

One way of comparing the various spheroidal popu- 
lations is to project them onto the CSph FP (Figure 
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Ej. We infer from the continuity of the various sam- 
ples in this plot that a single relationship may unite 
all systems from Es to CSphs. This underlying unify- 
ing relationship is more complex than a simple plane. 
The trend of decreasing A from ellipticals to BCGs to 
CSphs suggests that rather than three distinct planes, 
these spheroids may populate a slightly twisted 2-D 
manifold. To extend the range of spheroids, we also 
plot several other samples that contain lowe r a systems 
ijMatkovic fc Guzmanl2005HGeha et all2003|) . again cor- 
recting the pho tometry using ei t her m odel colors and the 
tabulations bv iFu kugita e t al.1 l)1995j) or measured col- 
ors (for the iGeha et al.l <|2003l) sample, we use the mean 
measured colors for dEs from lStiavelli et all (|2001)). Al- 
though the existence of a single, unifying manifold is less 
obvious for these systems because they deviate from the 
E-BCG-CSph trend, we show below that they can be 
included. 

We conclude this section by placing these systems in k- 
space (|Bender et al.lll99l . The CSph data appear as an 
extension of the elliptical and BCG data in all three pro- 
jections, again suggesting a physical connection among 
spheroids of different scales (Figure IHJl. In the k\ — k 3 
projection, which corresponds to the mass— log(M e /L e ) 
projection, the CSphs map onto the distribution of Es 
and extend that relationship up to the very large M/L 
ratios of clusters. The tightness of the trend in this panel 
suggests that much of the difference among spheroids 
reflects differences in the relative importance of dark 
matter wit hin r e . Previous investigators of the tilt 
of the FP llFaber et al.lll987t iPrueniel fc Simienl Il996t 
iTruiillo et alTf 2004) have considered the effect of system- 
atic changes in M e /L e with other galaxy properties and 
reached varying conclusions. The advantage provided by 
comparing our data to that of smaller spheroids is that 
the range of M e /L e is so large that some causes of M/L 
differences, such as stellar age and metallicity variations, 
are not physically plausible explanations for the entire 
relationship. Variations i n M/L with galaxy mass were 
previously identified (cf. IBurstein et aljfl99^) . but our 
analysis is different in its use of the p roperties of the 
CSph , rather than the member galaxies (Schacffcr et al. 
1993), and in the continuity and overlap among those 
samples, which enables cross-checking for systematic dif- 
ferences among samples. These advantages will enable 
us to identify a simple relationship between M e /L e and 
a. 

3.4. Beyond the Fundamental Plane 

The standard, back-of-the-envelope derivation of the 
FP begins with the virial theorem. Subsequent non- 
trivial simplifications, such as the use of a in the ki- 
netic energy term, the constancy of the numerical con- 
version factor between the measured line-of-sight veloc- 
ity dispersion and the 3-D velocity dispersion a for all 
spheroids, the use of r e in the potential energy term, the 
self-similarity of the potentials for all spheroids, as well 
as the assumption of sphericity, produce 

a 2 oc M e /r e , (2) 

where the proportionality constant depends on details 
of the structure that are assumed to be consistent 
among spheroids. The assumption that the structure of 
spheroids is scaleable and has no preferred size is referred 
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Fig. 7. — The projection of other spheroid populations onto the 
CSph FP. The extrapolation of the CSph FP to lower a is illus- 
trated with the solid line. Solid circle s represent our CSph mea - 
surements, open squares the BCGs of IQegerle fc Hoessell ft99ll) . 
op en triangles the Es o f p0rgcnscn et al (1996), crosses the Es 
ofjMatko«c ji£_ Cuzmanl (2205), and solid triangles the dEs of 
IGeha et al.l 120031 . There is no single FP for all spheroids. 
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Fig. 8. — Ellipticals, BCGs, and CSphs plotted in re-spacc. 
Solid circles represent our CSp h mea surements, open squares 
th e BCGs of IQegerle fc Hoesse] 119911) . o pen triangles the Es 
of IJprgensen et all 119961) . crosses th e Es oflMatk ovic fc Guzman! 
120051) and solid triangles the dEs of IGeha et alJ 12003). The line 
in the lower panel marks the boundary identified by Burstc in et all 
IT9971) . 



to as homology. By converting mass to (M/L)L, Eq. 2 
is rewritten as 

a 2 oc {M e /L e ){I e r 2 )/r e . (3) 

Taking the logarithm and rearranging terms, we arrive 
at 

log r e = 2 log a - log I e ~ \og(M e /L e ) + C, (4) 

which is the source of the expectation that the FP coef- 
ficients A and B should be 2 and 1, respectively (see Eq. 
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TABLE 3 

Direct Fundamental Plane Fits For Different 
Spheroids 



Components 
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1). The "tilt" in the FP, which corresponds to A and 
B not having the values this simple calculation suggests, 
is often attributed to the neglected M e /L e term, or to 
the breakdown of the various simplifying assumptions, 
principally homology. 

One difficulty in addressing whether Eq. 4 is an accu- 
rate description of spheroids is that there is no a priori 
expectation for how M e /L e varies with any of the other 
parameters. That dependence has generally been solved 
for after the fact by requiring it to reproduce the observed 
value s of A and B l|Faber et al.lll987tlPrugniel fc Simienl 
f996). This approach is underconstrained because a so- 
lution is guaranteed if a tilted plane originally fit the 
data. In other words, the reverse engineering of M e /L e 
is telling us either something about how it varies with 
galaxy properties or that the simple formulation of the 
FP in Eq. 4 has failed. External checks, such as a com- 
parison between the inferred M e /L e and the colors of the 
stellar population, suggest that there is at least some de- 
pendence between L and M e /L e . However, a wider range 
of data, including K band photometry, which is less sen- 
sitive to variations in stellar populations, demonstrates 
that only a small fraction of the overall tilt observed in 



the elli ptical FP ca n com e from stellar popula tion differ- 
ences (Eihreeral]|l99l ^uHIIoeranilOOl). Because 
the change in M e /L e across the range of normal elliptical 
galaxies is modest (typically a factor of a few at most), 
there are several different physical causes that could give 
rise to the FP tilt. Observing a sample that spans a much 
larger range of luminosities or velocity dispersions pro- 
vides a critical test of the inferences drawn from normal 
ellipticals alone f see iNapolitano et al.ll2005|) . 

On the basis of the trend in A values seen in Table 
Eland the tightness of the sequence between k1 and k3, 
we explore the connection between a and M e /L e . We 
estimate M e J L e using a dimen sional argument [M/L oc 
a 2 /(I e r e ); Burstein et al.lll997l) and plot the relationship 
between M e /L e and a in Figure El It is evident that the 
relationship between these two quantities is not a simple 
power law, as would be the case if the homology assump- 
tion were correct. The range in M e /L e over which the 
relationship is now evident rules out earlier interpreta- 
tions of apparently linear trends, seen only over a mod- 
est range of M e / L e , as arising from stellar population 
differences. While it is possible that stellar population 
differences account for some of the M e /L e variation, par- 
ticularly for systems with small M e / L e , they cannot ac- 
count for the factor of several hundred change in M e /L e 
seen here. The same can be said about suggested causes 
of homology breaking such as orbital anisotropy differ- 
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Fig. 9. — Me/Le versus a. M e /L e is calculated using K3 and 
is in solar units corresponding to the /—band (data other than 
ours is corrected to this band; see text). Because these values are 
calculated assuming homology, they are not internally consistent 
with our finding that M e /L e does not scale as a power law of 
a and may not represent in detail the true relationship between 
Me/Le and a. Nevertheless, the qualitative result of a non-linear 
relation between M e /L e and <r is robust. Solid c ircles represent our 
CSph measurements, open squa res the BCGs of lOegerle He Hoessej 
119911) , open triangles the Es of | j0rgensen et al 1199a) . crosses the 
Es of | Matkovic~S ^ Guzman (2005i), and solid triangles the dEs of 
IGeha et alj 120031) . 



ences or systematic differences in the surface brightness 
profiles. 

The variations seen here, if attributable primarily to 
one cause, appear to reflect the efficiency of concentrat- 
ing and converting baryons to stars within r e . The in- 
ability of a power law to describe the relationship be- 
tween M e /L e and a is clear from Figure El but it can 
also be easily demonstrated that no simple power law de- 
pendence of M e / L e on a can reproduce the variation in 
the FP A coefficient seen in §3.3. Consider M e f L e oc a a . 
Such a relationship changes Eq. 4 to 



logr e = (2 - a) logcr - log/ e + Ci. 



(5) 



While such a power law relationship could explain why 
FP fits do not find A = 2, it would not explain why A 
varies with a. Similarly, postulating power law depen- 
dences of M e /L e on I e and r e would change the numerical 
value of the coefficients in the FP, but not the functional 
form. Therefore, the description of M e /L e as a 2 /I e r e 
cannot account for the varying tilts of the FPs. 

Instead, from an inspection of Figure we suggest a 
new relationship: 



log M e /L e = {aloga-P) 2 



•7, 



(6) 



where a, (3, and 7 are free parameters. Such a relation 
modifies the original FP equation (Eq. 4) to 

log r e = (2 - a 2 log a + 2af3) log a - log I e + C 2 (7) 

and results in an effective A coefficient that is a function 
of a. One can imagine invoking other functional forms 
for M e /L e (a), but an offset parabola, a 2nd order poly- 
nomial, is the simplest mathematical form that can fit 
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TABLE 4 
Measurements Used in Fitting the 
Fundamental Manifold 3 
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the nonlinear relationship in Figure ED This form, which 
describes a twisting plane or a 2-D manifold, is not a 
unique choice, but we will show that the quality of the 
fit leaves little room for improvement with alternative 
fitting functions. 

3.5. Relationship between M e / L e and Velocity 
Dispersion 

We fit a parabola to the M e /L e data for dwarf ellipti- 
cals, ellipticals, BCGs, and CSphs in Figure It is im- 
portant to clarify two issues before proceeding. First, the 
M e /L e values that we fit are not derived from detailed 
dynamical modeling, but are rather those calculated from 
a dimensional argument. As such they cannot be the 
correct values in detail because, as we determined above, 
M e /L e cannot be described as a simple ratio of c, I e , and 
r e and produce the correct scaling of the A coefficient. 
Nevertheless, these M e /L e values should roughly track 
the true values because deviations from a single FP fit 
to all of these spheroids are not gross (rms scatter ~ 0.2, 
see below). Second, these values of M e /L e , as well as 
those necessary to fit Eq. 7, reflect the ratio within r e , 
not the global (total) M/L. 

In Figure ED we plot the same elliptical, BCG, and 
CSph M e /L e values versus a as in Figure |H1 but with 
several key additions. First, we have added the fit to 
the data in Figure El (solid line for the region constrained 
by data, dashed for the extrapolation). Second, we have 
adde d the data for dwarf ellipticals and dwarf spheroidals 
from | B ender et al. ( 1992) and for dwarf spheroidals from 
IMateol l|1997jh We have corr ected both data sets t o the 
/-band using the models oflFukugita et alJ 1)1995(1 and 
colors of dEs llStiavelli et al.ll20011) where necessary, and 
corrected the lBeridljr'l>rinT|l99*^ data using the cur- 
rently accepted distances to the Local Group dSphs and 
our adopted value of Ho for galaxies outside the Lo- 
cal Group. These resolved systems were not included 
in the fitting, or in Figure El because their properties 
(I e , r e and M e /L e ) are measured so differently than for 
the other systems that systematic errors in the compar- 
ison are likely. Even so, their properties are in concor- 
dance with the extrapolated parabolic fit to the proper- 
ties of the unresolved systems. While there is large scat- 



ter in the M e /L e measurements at the lowest c's, the 
data support our suggestion of a parabolic description 
for the relationship between logM e /L e and logc. The 
detailed fit cannot be used to quantitatively formulate 
the "fundamental manifold" because the underlying as- 
sumptions are not internally self-consistent. We exclude 
globular clusters from this discussion because they are 
not spheroids embedded in a dark m atter potenti a l well. 

We fit the data, fo r dE s iGeha et all 12003 
iMatkoyic fc Guzman] 120051) . Es l|j0rgensen et all fl996). 
BCGs i|Oegerle fc Hoessellll99lD and CSphs to a slightly 
rearranged version of Eq. 7, which we term the 
"fundamental manifold" in analogy to the FP, 

logr e = -a 2 log 2 (7+(2 + 2a/?)logCT + 5log/ e + C*2. (8) 
The resulting parameters are a 2 = 0.63 ± 0.02, (2 + 
2a/?) = 3.70±0.07,B = -0.705±0.004 and C = -2.75± 
0.08. If one uses only the ICS luminosity for the CSphs, 
then the resulting parameters are a 2 = 0.66 ± 0.02, 
(2 + 2a(3) = 3.81 ± 0.08, B = -0.676 ± 0.004 and 
C" = -2.90 ± 0.08. The comparable plot to Figure is 
shown in FigurclTTland the data used to obtain this fit are 
presented in Table 01 The rms scatter along the r e axis 
is now 0.099 for the combination of the Es, BCGs and 
CSphs (0.107 if the CSph is described by the ICS alone), 
which is not significantly worse than the 0.094 scatter of 
the elliptical FP alone. With the addition of the dwarf el- 
liptic al samples l|Geha et alJl2003l iMatkovic fc Guzmanl 
2005)), the scatter rises to 0.114. The combined sam- 
ple is highly weighted toward normal ellipticals, simply 
because they dominate the sample numerically, so it is 
important to examine the scatter carefully. For exam- 
ple, fitting a single FP to these same data results in a 
scatter of 0.129 between observed and fitted r e , which 
is not much worse than that derived from the manifold 
fit. However, the rms deviation for the CSphs about the 
the manifold is 0.095, while it is 0.207 for the single FP 
fit. At the other extr eme in the a range, the dEs from 
the lGeha et alJ l|2003|) sample have an rms deviation of 
0.165 about the manifold and 0.256 about the single FP. 
The manifold results in significant improvements at both 
ends of the a range. 6 

There are many reasons why our manifold fit must re- 
sult in a scatter larger than the intrinsic scatter: 

First, we combine data from various surveys carried 
out in different filters with different analysis techniques. 
Although we attempt to correct the colors to the same 
system, these corrections are based on mean observed 
or modeled colors for the population and mean fil- 
ter— (— detector transmission curves. 

Second, we do not correct for any wavelength depen- 
dence on r y., which has been shown to exist at least for 
ellipticals (Pa hre et alJll998|) . While this is a systematic 
error within a given sample, in the merger of samples 
presented here it would increase the calculated scatter. 

Third, we use the velocity dispersion of the 
cluster galaxies rather than the I CS. Theoreti- 
cal modeling l)Murante et all 120041: iWillman et alJ 

e A comparison of the scatter of CSphs, as described by either 
the ICS+BCG or ICS luminosity, about the fundamental manifold 
fits again favors the ICS+BCG description of the CSph (the scatter 
in logr e about the fundamental manifold fit, when allowing for a 
normalization error, is 0.095 for the ICS + BCG vs. 0.138 for the 
ICS, and 0.095 vs. 0.153, respectively, when not allowing for a 
normalization error). 
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Fig. 10. — An extended view of M e /L e versus a. To the data presented in Figured! we have added dEs and dSphs_fromJBenderet^y] 
11 993) f filled squa res) and dSphs (open circles) fromjMatco 1 1993). Five of the dSphs, at logo- ~ 1, appe ar in both thelBender et alj 11993) 
and lMated 119971) compilations, although the M e /L e values are calculated in independent ways. For the IMateol (119971) data, we adopt the 
lower M e /L e value if two are given. We also add a curve illustrating the bes t-fit parabola. We use a dashed line for logcr < 1.4 because 
the fit is unconstrained for these values. Neither the Bender ct al. (1992) nor Mateo (1997) data are used in the fitting due to the large, 
possibly systematic, uncertainties. The simple parabola fits the data for systems ranging from dSphs to rich galaxy clusters. 



l200i IS ommer-Larsen et alJ l2005j) and observations 
(|Zibetti ^1^012005^ conclude that the ICS is somewhat 
more centrally concentrated than the cluster galaxies, 
and hence could have som ewhat different kinematics. 
iSommer-Larsen et all l)2005j) suggest that at small radius 
the ICS could have about one-half the velocity dispersion 
of the galaxies. Indeed, a difference between the veloc- 
ity dispersion of the cluster g alaxies and ICS has been 
detecte d in the Coma clus ter l|Cerhard et al.l 120051) . al- 
though |Kelso^On <|2002|) find aics — * a at large radii 
in another cluster. If a is roughly a constant fraction 
or multiple of (JicSi then this difference would be ab- 
sorbed into the constant term of Eq. 8. If a is only 
crudely related to oics > then we would find a large scat- 
ter in our relationships (making those relationships more 
difficult to identify). If the relationship between a and 



a ics depends on other cluster parameters, such as the 
concentration, then this would introduce a tilt in our re- 
lationships and move the CSphs away from the BCG FP. 
The continuity between the CSph and BCG FPs argues 
against this being a large effect. In any case, we expect 
any problems in using the galaxy velocity dispersion as a 
proxy for the ICS velocity dispersion to be an additional 
source of noise. 

Fourth, we apply no correction for non-sphericity, or- 
bital anisotropy fcf. lCiotti et al.ll996l). or rotation in any 
of these systems fcf. iBusarello et alJlO Q?). While the ro- 
tation velocity of the ICS is presumably small in compar- 
ison to the velocity dispersion, the large ellipticity of the 
ICS components must at least necessitate a correction 
when converting from line-of-sight velocity dispersion to 
the velocity dispersion applicable for the FP. 
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Fig. 11. — Edge-on view of the 2-D manifold. Units given i n text. All lu minosities are corrected to the /-band (see text). Solid circle s 
represent our C Sph measurements, open squ ares the BCGs of lOegerle fc H pessell 119911) , open triangles the Es of j0recnscn et all 119961) . 
crosses the Es of lMatkovic Guzman (2005), and solid triangles the dEs of lOeha et alf 120031) . The line is the 1:1 line, and for this choice 
of axes it represents the best fit manifold. 



Fifth, we do not correct the measured velocity disper- 
sion to the velocity dispersio n within r e . As suggested by 
iSommer-Larsen et, alJ |2005) i an additional complication 
is that this correction may be different for the galaxies 
than for the ICS. 

Sixth, we ignore any direct dependence of M e / L e on 
I e or r e . To the degree that a can predict I e or r e , we 
will model variations in M e /L e with those parameters, 
but we know from the fact that spheroids do not fall on 
a fundamental line that a alone cannot reproduce the 
behavior in I e and r e exactly. In particular, we have not 
accounted for B =/= 1 in our scenario. 

Seventh, we adopt the parabolic description of the re- 
lationship between M e / L e and a as a mathematical con- 
venience. The correct functional form could be asym- 
metric, or have a flatter bottom than a parabola, or less 
curvature. All of these issues suggest that the intrinsic 



scatter about the true "fundamental manifold" could be 
smaller than that measured here. 

As with most explorations of the FP, one is left with 
intriguing relationships but no direct understanding of 
the physics that gives rise to those relationships. At the 
core of all of these fitting functions is the virial theorem, 
but it is the deviations from the simple expectations that 
can illuminate patterns in galaxy formation. In our case, 
we have two distinct results to understand: (1) why does 
M e /L e rise with increasing a for large a and (2) why 
does M e /L e rise with decreasing a for small cr? The bal- 
ance between these two, possibly distinct, physical causes 
will account for the existence of the minimum M e /L e 
and its particular value. In broad terms, and in differ- 
ent contexts, both of these trends have been noted be- 
fore (generally for global M/L), and various investigators 
have attempted to explain them (cf. iBabul fc Reeslll992t 
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[ Martini 119991 iBenson et al 112000) iMarinoni fc Hudson! 
120021 Ivan den Bosch et a,l] 1200ft IDekel fc Wod 1200311 

The behavior at low a has drawn more int erest because 
1) the observational effect s — high M e JL^ (lMateolll997l) 
and low effective yields ijGarnettl l2002t iTremonti et all 
12004^ — have been known for a longer time and 2) the 
formation hist ory of low mass halo s is critical in hierar- 
chical models llBabul fc Reesl ir9921 Two primary pro- 
cesses are cited as responsible for the comparatively low 
star formation efficiency in low mass systems: winds and 
UV photoionization. One expects that supernova winds 
will more easi l y remove material in lower mass systems 
l|Martinl Il999t iDekel fc Wool 12003ft and that the inter- 
gala ctic UV radiation fi eld will more easily ionize the 
gas llBabul fc ReesllT992l 

The behavi or on the largest sca l es has been explored 
less, althoueh lPadmanabhan et"aT] l)2004fl found an anal- 
ogous rise in M/L as a function of mass among Es, and 
various groups hav e recently begun modeling the for- 
matio n of the ICS llMurante et alJl2004t iWillman etTaD 
12001 ISommer-Larsen et alJ I2005D . Those simulations 
predict that the ICS forms from the mergers of groups 
early in the history of the cluster. Because the ICS 
baryons are already in the form of stars at that stage, 
they behave much like the dark matter during the merg- 
ers. Therefore, they do not concentrate toward the center 
of the potential well and do not form systems with low 
M e /L e . 

We conclude that the trend between M e /L e and a de- 
scribes how well baryons are able to dissipate and con- 
centrate toward the center of the potential. The various 
driving factors are the degrees to which 1) the baryons 
have already been converted to stars prior to the assem- 
bly of the system, 2) new gas can cool and form stars 
in the final system, and 3) winds expel gas and lower 
the binding energy of the central component. Theo- 
retical work is beginning to explore how mergers affect 
the relative distribution of dark and luminous matter 
(|Bovlan-Kolchin et alll2005|) . 

3.6. Interpreting the Minimum 

The choice of a parabola to describe the logM e /L e 
versus loger relation implies a well-defined global mini- 
mum in the relationship; however, we have stressed be- 
fore that this particular functional form is somewhat ar- 
bitrary. Even so, the data, particularly for the dSphs, 
do suggest an upturn in M e /L e at low a, although a flat 
trough is not ruled out. If we accept the assumption of 
a parabolic relationship between M e /L e and a, our best 
fit manifold from Eq. 7 implies a — 0.79 ± 0.03 and 
= 1.07 ± 0.05 and that the minimum M e /L e occurs 
in systems with a = 22.9 ± 4.7 km sec -1 . The degree 
to which these values are given physical significance is in 
large part determined by the validity of the mathematical 
form we have chosen, how well the various samples have 
been placed on the same photometric and kinematic sys- 
tems, and the sampling of the data through the region 
of the minimum. Unfortunately, the bulk of our data 
are at larger a than the minimum and so our fit is quite 
sensitive to the limited data available near the minimum. 

In contrast to what we did above, the turnover in the 
M e /L e resulting from dimensional arguments (§3.4) oc- 
curs at a ~ 62 km sec -1 (Eigure lTU)) . While we have ar- 
gued that the evaluation of M e /L e using homology and 



dimensional arguments is not internally self-consistent 
with our modeling of the fundamental manifold (because 
of the higher order dependence of log M e /L e on logcr 
that we impose), these estimated values of M e /L e do 
reflect something about the relationship between a, r e , 
and I e . The reason the estimates do not agree in de- 
tail is presumably because the correct coefficients are 
not available in the dimensional argument and because 
the homology assumption breaks down. Which turnover 
is the more relevant in assessing the relative effects of 
winds, or photoionization, or structural changes due to 
the degree of dissipation is unclear. We conclude that 
while a change in behavior does occur, the exact value 
of a ranges between at least 20 to 60 km sec -1 . The 
corresponding mass scale is much lower than that sug- 
gested by studies of the global M/L that also find a min- 
imum Ml L at a specific m a ss scale (IBenson et al. 1120001 
IMarinoni fc Hudson! 120021: Ivan den Bosch et alJ 12003^ 
but the results are not necessarily in conflict. Observ- 
ing differences in the behavior of M/L within different 
radii, in this case within r e versus r2oo i could potentially 
be used to break the degeneracy in the analysis between 
the dissipation of baryons and the efficiency of converting 
them to luminous matter. 

3.7. Caveats 

As with all analyses of the FP and its cousins, the 
description of spheroids as fully virialized, homologous 
(or at least weakly homologous) dynamical systems is 
oversimplistic. Our hope is that a minimalist description 
will embody the principal elements of spheroids, and the 
mere existence of the FP supports this approach. Never- 
theless, to ignore the richness of the process of spheroid 
formation could lead to an overinterpretation of fortu- 
itous coincidences. Even in the case of low scatter in the 
FP (~0.1), the scatter between observed and model r e 
is still 26%, which, while small considering the range of 
r e values, is not negligible. For the analysis and discus- 
sion we have presented, we remain concerned about the 
following: 

We have ignored possible variations of M e /L e that de- 
pend on r e and I e . The standard approach attributes the 
tilt in the elliptical FP as being partly due to a and partly 
to the other terms in Eq. 4. We justify placing the entire 
burden on the form of the a term on the basis of the large 
variation in A among FPs and the success of the model. 
However, given the various correlations between r e , I ei 
and a it is possible to slough off some of that variation 
to the other parameters. In a complete treatment, one 
would aim to explain the behavior of the B coefficient in 
Eq. 8 as well as that of A. However, examining the resid- 
uals from the best fit manifold versus I e , we find only a 
marginal indication of systematic behavior, and hence 
expect a minimal reduction in the scatter with the addi- 
tion of a 7 e -dependent term. Instead, we choose to limit 
the number of new parameters introduced at this stage, 
but, as samples increase, particularly for extremely low 
and high a spheroids, we anticipate more complex and 
complete modeling. 

We have ignored other sources of homology breaking in 
our analysis. For example, systematic changes in orbital 
structure or in the dark matter potential as a function of 
spheroid scale will also tilt the FP. We argue that such 
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changes cannot account for the large (factor of ~ 100) 
difference in our calculated M e / L e across the entire range 
of spheroids, but they are presumably there at some level. 
To the degree that they vary as a low-order function of 
cr, they are incorporated into Eq. 8. 

We have focused on using the combination of BCG + 
ICS luminosity to describe the CSphs. In no other type 
of spheroid did we explicitly combine two components. 
Are CSphs truly described best by the combination, or 
are we being misled by small number statistics? Or do 
other systems, such as normal ellipticals, actually have a 
low surface brightness outer component, analogous to the 
ICS, that has not been observed? Regardless of whether 
we study the ICS alone, or in combination with the BCG, 
the existence of a fundamental manifold and the system- 
atic behavior of M e /L e with cr remain (Figure E} - 

The trends in A and in M e / L e are dominated by inter- 
sample differences rather than by intra-sample differ- 
ences. Because the samples come from disparate analy- 
ses, we are concerned that some of the difference among 
them is due to systematic errors rather than physical 
causes. Because there is significant overlap between sam- 
ples, we conclude that this is not a serious problem, but 
it is a fact that each of the principal samples alone is ac- 
curately described by a F P relation. Marginal e vidence 
for curvature was noted by Jorgc nsen et all l)1996j) in their 
sample, although, Bernardi et alJ l )2003|) accounted for all 
the apparent curvature in their sample through selection 
effects. 

Because most of the available data is on the high a side 
of the M e /L e minimum, we have little constraint on the 
functional form of the M e /L e — a relation. For exam- 
ple, an equally reasonable attempt at unification could 
come by modeling the M e /L e variation as a function of 
enclosed mass, r e a 2 . This approach leads to an equation 
that is somewhat harder to interpret because of cross- 
terms, but that has the same number of free parame- 
ters. Even if the functional form, whatever it may be, is 
steeply rising on both sides (as reinforced by the dSph 
data), there is no reason to believe that the rise should 
be symmetric. Once one allows for asymmetry, the po- 
sition of the minimum becomes quite uncertain and the 
need for data at logo - < 1.5 acute. 

The interpretation of r e depends on the surface bright- 
ness profile. There is still some uncertainty in the liter- 
ature regarding the relative merits of Sersic versus de- 
Vaucouleur profiles across all spheroid classes. The lat- 
ter type of profile is a special case of the former, so the 
more general statement that Sersic profiles fit spheroids 
is less controversial. A concern here is that there is a 
systematic trend in profile shape with a that produces 
correspond ing systematic error s in r e and I e , and hence 
in the FP. iBertin et al.l l)2002|) have explored this topic 
and concluded that "weak" homology is satisfied, but 
the exact effect on a n ana lysis such as ours is unclear. 
iGraham &: Guzman! ((2003) find a continuous relation- 
ship in structural parameters of dEs to Es when they 
fit Sersic profiles rather than r 1 / 4 profiles. Which pro- 
file is more appropriate across the full range of spheroids 
and how to compare among objects, different studies, 
and different methods are open questions. However, pro- 
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Fig. 12. — Same as Figures fTTTl and ITTI except the CSph lumi- 
nosity is taken to be only that of the ICS, rather than that of the 
ICS+BCG. In the left pane l, th e solid line is the best fit parabola 
to the data shown in Figure HHI and is included here as a reference. 
In the right panel, Y = 3.81 log 2 cr - 0.66 log a - 0.68 log I e - 2.90, 
an d rep resents the best fit surface. Symbols are as described in Fig- 
ure [IT]] Whether the CSph is described by the ICS or ICS+BCG 
does not affect either the relationship between M e /L e and a or the 
nature of the fundamental manifold. 



file differences seem unlikely to result in the large FP 
differences observed among the full range of spheroids. 

4. CONCLUSIONS 

We demonstrate that spheroids ranging from dEs to 
the intracluster stellar component of galaxy clusters, 
which we designate the cluster spheroid (CSph), lie on a 
curved surface, a 2-D manifold, in (cr, r e , I e ) — space. Pre- 
vious studies have generally examined limited regions of 
parameter space, over which the manifold can be ade- 
quately described by a plane. Our principal findings are: 

• The addition of cluster velocity dispersions, as mea- 
sured from the cluster galaxies, decreases the scatter of 
the r e -I e relation for the cluster spheroids, and there- 
fore signifies the existence of a "fundamental plane" for 
clusters. 

• Combining the luminosity of the brightest cluster 
galaxy (BCG) and the intracluster stars (ICS) leads to 
significantly lower scatter in the CSph FP relationship. 
We argue that this effect is real, implying a connection 
between the evolution of the BCG and ICS, rather than 
a poor decomposition of the two components. Our re- 
sults do not change qualitatively if we use only the ICS 
for the CSph, but the scatter in various relationships in- 
creases significantly. Unless otherwise noted, we refer 
to the combination of BCG and ICS as the CSph, but 
in the more important cases present results for both de- 
scriptions of the CSph. 

• The scatter of the CSph FP (0.074) is as small as that 
observed for elliptical galaxies, but the orientation of the 
plane is different. 

• There is a systematic decline in the coefficient of the 
logo - term in the equation of the FP as one progresses 
from systems with smaller to larger cr. This trend sug- 
gests that a 2-D manifold, which over the limited range 



Fundamental Manifold of Spheroids 



15 



of a probed by any one spheroid population is indistin- 
guishable from a plane, roughly forms a set of twisted 
planes that fit spheroids from low mass ellipticals to the 
CSphs of the most massive clusters. 

• The non-linear relation between log Me/ L e and a for 
spheroids implies a functional relation between these 
quantities that is more complex than a power law, and so 
requires the breaking of the homology assumption. We 
adopt log M e /L e = (a log <7 — (3) 2 + 7 as a model of this 
relationship because it reproduces the systematic decline 
in the coefficient of the log a term in the FP equation and 
captures the curvature seen in the logM e /L e versus logo- 
relationship for spheroids. Of particular importance is 
that M e /L e represents the mass to light ratio within r e 
and may be decoupled from the global M/L. The mag- 
nitude of the change in M e /L e across the full range of 
spheroids excludes changes in stellar populations as a 
possible explanation for the global trend. Instead, these 
changes must be driven primarily by variations in the 
relative amounts of luminous and dark matter within r e , 
although some modest dep endence of M P / L e with s tellar 
populations must exist fsee lCappellari etafl l)2005|) for a 
recent exploration of this topic for E's). 

• The resulting 2-D manifold has rms residuals along 
the r e axis of 0.099 (for Es, BCGs, and CSphs), com- 
pared to 0.094 for the FP plane fit to ellipticals alone 
and to 0.089 for BCGs alone. If we extend the mani- 
fold to dEs, the scatter increases slightly (0.114). The 
increased scatter relative to the individual FP's of each 
spheroid population is insignificant given that the mani- 
fold fit is done with respect to five different samples that 
were observed in different photometric bands by different 
investigators. Using the best fit 2-D manifold to define 
the coefficients of the relationship between logM e /L e 
and logcr, we find a minimum M e /L e for galaxies with 
a = 22.9 ± 4.7 km s _1 . We qualitatively agree with a 
range of theoretical and observational studies that con- 
clude that the re exists a galaxy mas s scale for which M/L 
is minimized l|B enson et al. I l200Ct iMarinoni fc Hudson! 
l2002Hvan den Bosch et al.ll2003ft . although our value for 
that mass scale is significantly lower. Part of the dis- 
crepancy may arise from the measurement of this phe- 
nomenon on global scales, r2oo> versus within r e . 

Our new description that places the family of spheroids 



on a 2-D manifold in (cr, r e , I e ) — space is remarkable in 
that it unifies spheroids that span a factor of 100 in a 
and 1000 in r e . This range of systems must have dif- 
ferent histories both in terms of when they were assem- 
bled and when their stars formed. Even so, the global 
properties of these systems end up on the manifold with 
intrinsic scatter that must be less than 0.114 (30% in 
r e ). We are able to place all dark matter dominated 
spheroids on a single relation by accounting for the ob- 
served non-linear relationship between log(M e /L e ) and 
log a\ This relationship, which deviates from simple 
virial theorem and homology expectations, points to a 
more specific phrasing of the question of how galaxies 
and clusters form: what is the origin of the relationship 
between M e /L e and potential well depth? The relative 
importance of dissipation, cooling, star formation, feed- 
back, and dynamical relaxation during spheroid forma- 
tion must vary with spheroid scale in such a manner as 
to reproduce the continuous trend observed in M e /L e 
(Figure IUH , a trend responsible for the "fundamental 
manifold" (Figure fTTf) . Theoretical explorations, such as 
those focusing on the interplay between different feed- 
back me chanisms and the growth of galaxies (for exam- 
ples, s eelPadman abhan et all2004l iBovlan-Kolchin et al 
120051: iDekel k, Brinboiml |2005|) . may ultimately reveal 
the nature of the observed trend in M e /L e and thereby 
unify the formation processes of spheroids. 
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